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Abstract: Axial charge imbalance is an essential ingredient in novel effects associated 
with chiral anomaly such as chiral magnetic effects (CME). In a non-Abelian plasma with 
chiral fermions, local axial charge can be generated a) by topological fluctuations which 
would create domains with non-zero winding number b) by conventional non-topological 
thermal fluctuations. We provide a holographic evaluations of medium’s response to dy¬ 
namically generated axial charge density in hydrodynamic limit and examine if medium’s 
response depends on the microscopic origins of axial charge imbalance. We show a local 
domain with non-zero winding number would induce a non-dissipative axial current due to 
chiral anomaly. We illustrate holographically that a local axial charge imbalance would be 
damped out with the damping rate related to Chern-Simon diffusive constant. By comput¬ 
ing chiral magnetic current in the presence of dynamically generated axial charge density, 
we found that the ratio of CME current over the axial charge density is independent of the 
origin of axial charge imbalance in low frequency and momentum limit. Einally, a stochas¬ 
tic hydrodynamic equation of the axial charge is formulated by including both types of 
fluctuations. 
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1. Introduction 


The parity-odd response of a medium with chiral fermions and its deep relationship to 
topology and quantum anomalies have attracted significant interest. One such effect under 
extensive study is the chiral magnetic effect (CME) [1-4] , which is the appearance of a 
vector current along the direction of an external magnetic field in the presence of axial 
charge imbalance (see Refs. [5,6] for a recent review). The CME has been demonstrated in 
various theoretical frameworks, such as in hydrodynamics [7-13], kinetic theories [14-19], 
perturbative theories [4,20,21], effective theories [22-26] and in the AdS/CET correspon¬ 
dence [27-32]. A closely related effect is the chiral vortical effects (CVE), which is the 
appearance of a current along the direction of vortivity. CVE and its relation to mixed 
anomalies has been studied in [33-37]. Interesting properties of chiral media have been 
discussed in [38-40]. 

Those anomalous effects are not only theoretically well-motivated, but also phenomeno¬ 
logically important. In a heavy ion collision, a very strong magnetic field, on the order 
of eB ~ m^, is created from the incoming nuclei that are positively charged and move 
at nearly the speed of light. Therefore, CME will convert axial charge fluctuations gen¬ 
erated in heavy-ion collisions into (vector) charge-dependent correlation which could be 
potentially detected by experimental observables. Recently, there have been significant 
experimental efforts in searching for CME and other anomalous transport effects (see [41] 
for a review) in heavy ion collision experiments [42-45]. 

One essential ingredient in those anomaly-related effects is the presence of axial charge 
imbalance. For example, in terms of chiral charge imbalance parametrized by the axial 
chemical potential CME can be expressed as; 

= CAfiAeB, Ca = ^. ( 1 . 1 ) 

Previously, most studies were based on introducing axial charge asymmetry by hand, after 
which the response of the medium to a magnetic field is investigated (see Ref. [46,47] for 
exceptional cases). However, axial charge density nA{t, x) is a local and dynamical quantity 
depending on space and time. The medium’s response to time-dependent, in-homogeneous 
axial charge density and the connections of this response to anomaly have been rarely 
studied before. One motivation of this paper is to fill this gap. 

A distinctive feature of local axial charge density (in contrast to vector charge density) 
is that there are microscopically two different mechanisms for local generation of axial 
charge imbalance. The first one is by topological fluctuations of gluonic fields which would 
create domains with non-zero winding number. The resulting topological charge will in 
turn convert into axial charge density via anomaly relation: 

df.fA = -‘^q: (1.2) 

where q is topological charge density (Pontryagin density) and denotes field strength of 
gluonic fields. Such topological fluctuations would in general create both global axial charge 
and local axial charge imbalance. The second mechanism is through thermal fluctuations. 
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In the absence of topological fluctuations, while the whole system is a grand canonical 
ensemble (of axial charge), each fluid cell of that system can be considered as a canonical 
ensemble (of axial charge). Therefore axial charge inside a fluid cell still fluctuates due 
to thermal fluctuations. For a non-Abelian plasma with chiral fermions, both mechanisms 
would contribute to local axial charge density fluctuations. Would response of the medium 
depend on the way that axial charge density is generated? 

In this paper, we will consider a de-confined non-Abelian plasma with chiral fermions 
at finite temperature. We will begin by studying the medium’s response to the interplay 
between local axial charge density nA{t,x) and jAit,x), q{t,x) in long time and long 
wavelength limit (i.e. in hydrodynamic regime). In particular, we will examine if the 
relation among those one point functions depends on the microscopic origin of local axial 
charge density. As a basis for this study, we will work in a top-down holographic model, 
namely Sakai-Sugimoto model [48,49]. The Sakai-Sugimoto model is considered to be 
close to the large W QCD with massless chiral quarks in quenched approximation. It has 
been widely applied to study anomaly-related effects (e.g. [27,50,51] ). To model gluonic 
fluctuations and implement anomaly relation (1.2), we will consider the dynamics of Cy 
Ramond-Ramond field and its Wess-Zumino coupling to flavor sector. As we are working 
in the long time, long wave length limit, the results presented in this paper are analytic. 

We found axial current jA{t,x) can be created by in-homogeneity of topological do¬ 
mains. It is well known that in-homogeneity of axial charge densities leads to diffusion: 
jA{t,x) = —D'VnAit,x) where D is the conventional diffusive constant. However, if the 
axial charge density results from a local topological domain which will be represented by 
an effective field “0(t, *)” in this paper , such topological domain will also induce an axial 
charge density and non-dissipative axial current in addition to diffusive current: 


r pQ 

UA = x) 


jT = «csV0. 


(1.3) 


Here, Tcs and kcs is related to the behavior of retarded Green’s function G|j(t, 
([(?(t, £c), g(0,0)]) in hydrodynamic regime: 


X 




.Tcs ,2 


(1.4) 


and T denotes temperature. It is worthy noting that kqs term in (O) is opposite to the 
direction of diffusive current and non-dissipative. One way to understand current (O) is 
that a topological domain also carries kinetic energy which would be transfered to chiral 


fermions via anomaly relation (1.2). In Ref. [52] by us, we have derived (1.3) based on a 
generic setting and presented a brief verification of (|1.3|) in the Sakai-Sugimoto model. We 
provide more details on this calculation in Sec. and elucidate how axial current (1.3) is 
generated from gravity side of the duality. 

Our calculation in Sec. |3.2| confirms that a local chiral charge imbalance nA{t,x) will 
induced a non-zero q{t,x) and they are related by 


q{t,x) = 


nA{t,x) 


2t, 


sph 


(1.5) 
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Tgph here can be interpreted as the axial charge damping time. Indeed, substituting (|L^) 
into (|1.2D , one has: = —UA/Tsph- Eqs. (|1.5D hence implies that a non-zero axial 

charge density will eventually be damped out by inducing a non-zero q. Furthermore, we 
verify by holographic computations that Tgph is related to Chern-Simon diffusive rate Fes 
and susceptibility y: 


xT 

’■•pi' - 2rcs ■ 


( 1 . 6 ) 


Previously, relation (1^) has been derived based on the standard fluctuations-dissipation 
argument [53] (see also Sec. ^). Very recently, Tgph has also been computed numerically in 
a bottom-up holographic model [54]. To best of our knowledge, current work is the first 
direct verification of relation (|1.6| ) in strong coupling regime. 

As we find that the axial current in response to axial charge density depend on how 
such axial charge imbalance is generated, it is natural to ask if chiral magnetic current (El) 
also depends on the origin of axial charge imbalance. To be quantitative, we consider the 
ratio between CME current and axial charge density in low frequency, small momentum 
limit in the presence of constant magnetic field: 


(Xdyn) — 


lim 

>-0 


nA(w,/c)J 




j^^^{uj,k) 

CacB 


(1.7) 


For a system with constant axial charge density n^, the ratio /{CeB) equals to 

axial chemical potential ^a due to (^). However, if ua is space-time dependent, the 
dehnition of axial chemical potential fiA is ambiguous. If one takes the ratio /{CeB) 

as the generalized definition of axial chemical potential, the ratio nAi^^,k)/fj,Aii-o,k) can 
be interpreted as susceptibility. For this reason, we will call Xdyn the “dynamical axial 
susceptibility”. 

We would like to emphasis that Xdyn (1-7) is conceptually different from chiral mag¬ 
netic conductivity [55] which is the proportionality coefficient of CME current to the time- 
dependent magnetic held for a medium with homogeneous, time-independent axial chemical 
potential. In ( |l.7| ) however, magnetic held is constant while ua is space-time dependent. It 
is worthy noting that in realistic situations such as quark-gluon plasma (QGP) created in 
heavy-ion collisions, Xdyn would be a relevant measure of CME as in those situations, the 
axial charge density is always generated dynamically. Previously, chiral magnetic conduc¬ 
tivity has been calculated for plasma in equilibrium at both weak coupling limit [20,21,55] 
and strong coupling [50,56,57] and for plasma out-of-equilibrium [58]. However, we are 
not aware any existing literature discussing the “dynamic axial susceptibility” Xdyn and its 
universality. 

We have computed Xdyn with nA{uj,k) generated by topological and thermal (non- 
topological) fluctuations in Sakai-Sugimoto model. We found such ratio Xdyn (0) is in¬ 
dependent of the origin of axial charge imbalance and equals to static susceptibility x- 
Moreover, we derive a simple analytic expression ( 4.15| ) relating the (integration of) gravity 
metric to Xdyn which applies to a large class of holographic model. From such expression, 
we obtain a condition on the universality of Xdyn- 
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Having established the fact that axial charge generate by both topological fluctuations 
and thermal fluctuations would contribution to CME current, it is then important to incor¬ 
porate both fluctuations in the framework of stochastic hydrodynamics. Recently, there are 
encouraging progress on applying anomalous hydrodynamics to simulate charge seperation 
effects [59, 60] and chiral magnetic wave effects [61, 62] in heavy-ion collisions. In those 
studies, axial charge density enters as the initial conditions while the fluctuations of axial 
charge density during hydrodynamic evolution have been neglected. Motivated by hndings 
in this paper, we formulate a stochastic hydrodynamic equation of axial charge density 
in Sec. Such hydrodynamic equation includes stochastic noise from both topological 
fluctuations and thermal fluctuations. While it is a direct generalization of the general 
framework [63-65], to our knowledge, stochastic equation (^) is new in literature. We 
hope our theory would be applied to simulate phenomenology of anomalous transport in 
the future. 

The paper is organized as follows. We will begin with a brief review of pertinent 
ingredients of Sakai-Sugimoto model and realization of anomaly relation ( 1 . 2 ) in Sec. |^. 
Sec. 1^ is devoted to studying medium’s response to axial charge density. The computation 
of Xdyn (0) is presented in Sec. The stochastic hydrodynamic equation for axial charge 
is formulated in Sec. We conclude in Sec. 


2. Sakai-Sugimoto model and chiral anomaly 


2.1 Set-up of the model and realization of anomaly 

In this paper, we will work in Sakai-Sugimoto model [48,49]. In this model, the de-conhned 
phase is given by DA black-brane metric, which is a warped product of a 5d black hole and 
X [66,67]. The DA brane background is given by [ 66 ]: 


ds^ = 


/r/\ 3/2 

f j f {U)dt^ dx^ dxl) -|- 




dt/2 \ 

mJ ’ 


( 2 . 1 ) 


z^RR 27riVce4 

D - 



= TrgsNJs , 


fiU) = 1 - 



( 2 . 2 ) 

Here X 4 is the coordinates of and €4 is the volume form of the four sphere S 4 . In 
addition, V 4 = Stt^/S is the volume of S 4 and gs,ls are string coupling and string length 
respectively. The location of the horizon Ut is related to the inverse temperature: 


47rR3/2 ^ I 


The periodicity of X 4 is given by 


6t = 2ttR4 


2 -k 


(2.3) 


(2.4) 
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The background (^) is stable for T > l/(27ri?4) = Mkk/‘^'^ [67]. Finally, ’t Hooft 
coupling A is given by: 


X = N, 


{iTrygsls 


= iTrNcgshMxK ■ 


(2.5) 


27ri24 

To model gluonic fluctuations, we will consider the dynamics of Cj Ramond-Ramond 
form. The kinetic energy of C 7 are given by 

{2Trlsf 


Srr = -- 


dvr 


[ dC7Ai*dC7). 
JlO 


( 2 . 6 ) 


In Sakai-Sugimoto model, right and left handed quarks are introduced by Nf D 8 branes 
and Nf D 8 branes [48,49]. Right-handed (left-handed) U{1) gauge filed Ar (Al) lives on 
D 8 {D 8 ) and is dual to right-handed (left-handed) current (J^) on the boundary. The 
D 8 /D 8 branes are separated along the X 4 direction, with D 8 branes located at X 4 = 0 and 
D 8 branes located at X 4 = 7 ri? 4 . In this work will consider Nf = 1 though generalization 
to the case of multi-flavors is straightforward. 

The action of bulk gauge held or its held strength is given by the summation 
of Dirac-Born-Infeld (DBI) term 


nR,L _ 1 

“ (27r)8/9 


j d^xe "^^J-det {gMN + (27ra') , (2.7) 


and Wess-Zumino (WZ) term which couples to Ramond-Ramond form; 


qR,L/ _ I 

^ 


/ 

JE9 


A tre 27 t 


( 2 . 8 ) 


We normalize the RR forms C as in [48] and use hermitian worldvolume gauge held A. The 
DBI action is identical for An and A^. In WZ term, plus/minus sign is corresponding to 
Ar!A n (i.e. F^j^/F^j^) respectively. Axial gauge held and vector gauge held are related 
to right-handed and left-handed gauge held by: 


A = 


AR-A^ 


V = 


AR 


The total action we will study then becomes; 

S = Srr -I- 5dbI + -^wz 


(2.9) 


( 2 . 10 ) 


It is instructive to show how axial anomaly relation (1.2) is realized in the current 
holographic model. Following holographic dictionary, the axial current is given by the 
variation of holographic on-shell action Shoio with respect the boundary value of = 
Afi{U —)• 00 ) and q is given by the variation of 5hoio with respect to 9: 

dSYiolo 


■ U '^‘S'holo 

5a,, 


Q = 


59 


( 2 . 11 ) 


Here 9 is determined by the holonomy of Ci on the compactihed X 4 direction [ 68 ]: 

9{t,x) = lim [ dx 4 {c[^^). (2-12) 

[/—>00 / 
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One may note that the normalization in (2 .12 ) is consistent with the one found by consid¬ 
ering action of multiple probe color branes. 

Cl is related to Cj and held [69]: 


dCi = {2'kIs)'° * dC^ — A { 6 {x 4 — 'KR/i)dx 4 ) + /\ {5{x/i)dx4) . 


(2.13) 


It is clear from ( 2.15 ) that Cj is invariant under axial gauge transformation: 

Sp^A^ = dA, dpiA^ = —dA, dp^Ci = —K5{x/i)dxi — K5{xi — 'KRi)dx/i , (2.14) 

where A is an arbitrary scalar function. At boundary and after integrating out Xi depen¬ 
dence, the axial gauge transformation is reduced to 


x) = x ), dpiOit, x) = -2A(t, x) 


(2.15) 


Since the action (ED expressed in terms of C^ and field strength Fr/Fl are manifestly 
invariant under axial gauge transformation ( p.l4 ), 5hoio would also be invariant under 
( 2.15 ). We therefore have that for an infinitesimal transformation 6A: 

‘^A'S'jiolo — j d x 


.d^{6A{t,x)) + ^^^{-26A{t,x)) 


5a^(t,x) ^ 69{t,x) 

J d‘^x x) {d^5A{t, x)) - 2q{t, x)6A{t, x)] = - J d^x x) + 2q{t, £c)] 6A{t, x) = 0. 


(2.16) 


The anomaly relation (O) then follows from the requirement that ( ^.16 ) holds for arbitrary 
(5A, This is the holographic realization of axial anomaly in the current model. Realization 
of axial anomaly in general Dp/Dq brane can be found in [70]. 


2.2 Fluctuations of bulk fields 

We wish to study medium’s response to local axial charge imbalance. To model that process 
in holography, we need to introduce sources on the boundary. Those sources will excite 
bulk fields which in turn would generate one point functions such as axial current and 
q on the boundary. As we discussed in the introduction, local axial charge imbalance can 
be generated by a gluonic configuration with non-zero winding number and by thermal 
fluctuations. Correspondingly, we will create axial charge imbalance by putting a non-zero 
9{t, x) and by putting a non-zero axial gauge field a^(t, x) on the boundary. In this work, 
we will restrict ourselves to the longitudinal fluctuations, i.e., the source on the boundary 
are 9{t,x),at{t,x),ax{t,x) where x corresponds to the direction of non-vanishing current 
and sources will only depend on t, x. Consequently, non-zero bulk fields are At, Ax, Ajj, Ci 
and they would only depend on t, x, U^. For later convenience, we introduce a dimensionless 
radial coordinate: 

u=UIUt. (2.17) 

^Note that C 7 can also depend on X4. This can happen when we consider backreaction of D 8 brane. We 
will restric ourselves to lowest mode on 5"^ with no X 4 dependence in this work. 
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We will also rescale all other dimensionful quantities by 


Ut dvr 


A-^h. 

3 


(2.18) 


where we have used 

We take following ansartz for Cy: 


C^ = BM{t, X] u) dx^ AdcJi 2 Ae 4 = x; u)dt + Bx{t, x; u)dx + Bu{t, x; u)du] A(icJi 2 Ae 4 . 

(2.19) 

where dai 2 = dxif\dx 2 - Here, M, N run over t, x, u. Adopting the convention eta;ix2X3X46»i6»26»36»4i 
1 with 6i(i = 1 , 2 ,3,4) being angular variables on S^, we obtain the explicit expressions for 
dC^ and *dCT. 

dCj = [Gtxit, x\ u)dt Adx + Gutit, x\ u)du f\dt + Guxit, x\ u)du A dx\ A dai 2 A 64 , (2.20a) 


*dCf = —u ^ [G^“(t, X] u)dt + G^^{t, x; u)dx + G^^{t, x; u)du\ A (ix 4 , 


( 2 . 20 b) 


where Gmn = QmBj^ — OnBm is the field strength of B. Indices are raised/lowered by 
metric (P), i.e. : 


gtt{u) = -u^^^f{u ), gxxiu) = , guu{u) = u ^ 


( 2 . 21 ) 


and f{u) = 1 — u Action ( 2.1C1| ), after performing trivial integration of X4, can then be 
written in a compact form: 

S = J d^xdu (^-^N.^F^^Fmn -\v^G^^Gmn - Kc^^^BlFmn^ , ( 2 . 22 ) 


where Fmn = Q^Am — OmAn is the field strength of A. Three terms in ( 2 . 22 ) are corre¬ 
sponding to Sdbi, Srr, Swz respectively and the sign convention for Levi-Civita symbol 
is = 1. In 


= CgU-^. = G^u^'^ . (2.23) 

The dimensionless parameter Cg , can be expressed in terms of Mkk and T as 

729ttK^MIj^ „ 2 AT 


c,= 


= 


(2.24) 


4A3r2 ’ ' 277rKMKK 

Varying (p.22) with respect to A and B, we obtain equation of motion for 


duiV^G^^) = Ke^^^FQM 


(2.25a) 


Nc {V^dMF^^) = Ke^^^GQM 


(2.25b) 
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2.3 The prescription for computing one point function 

In this section, we will derive the explicit holographic prescription for computing one 
point function nA,jAiQ from definition ( 2.11 ). For this purpose, we need to first obtain 
holographic on-shell action and then perform the variation with respect to the sources, 
9,at,ax- It is therefore more natural to rewrite the holographic action in terms of A and 
Cl field instead of Cy. The reason is that Ci is directly related to the boundary source 6 
as seen in Eq.( p.l2 ). From (|2.I3|) and ( p.l9|) , we have explicitly: 

= 2K [OlM + 2Al] . (2.26) 

where M is the X 4 component of Ci, integrated over the X 4 circle: 


M= dx^C 


(4) 


(2.27) 


The boundary value of M is 9. Note that and {d^M + 2 An) are invariant under 

both flavor and Cy gauge transformations (2.14). Using ( 2.26| ), ( |2.25b| ) can be written in 
terms of M as 

NcBm {d^M + 2A^) 

where we have defined: 

^ ^_ {~9tt9xx9uu) 




(2.28) 

(2.29) 


Moreover, the Bianchi identity of Gmn reads 

e^^^dLGMN = 0, 


together with relation (2.26) gives 


On 


( 271 ^ + d^M) 


= 0 . 


(2.30) 


(2.31) 


It is easy to see the action in terms of the M and Am fields, which would lead to equation 
of motion ( 2.28 ) and ( 2.31 ) is 


S= d'^xdu 


- \N^y/^F^^FMN - ^V^idqM + 2 Aq) {d^M + 2^^) 


(2.32) 


as it is analyzed in [69], [70]. The action has the same form as Eq. (9) in [71] from a 
bottom-up model. However our action differs in the Nc dependence of different terms in 
( 2.32 ), which is absent in [71]. 

In order to compute the one point functions of q,nA and Ja we use the action ( 2.32| ) 
instead of ( 2 . 22 | ), since it is expressed in terms of the bulk fields which are directly related 
to the sources of the boundary operators. To obtain the on-shell action, we do variation of 
(|0^ : 

6 S= [ d^xdu{dM INc^F^F^^ - 2K‘^^/^{d^M + 2A^)] 6 An + K^y/^dN {d^M + 2A^) 5M} 


+ / d‘^x 


Ncyf^F^^dAN - + 2A^)6M 


(2.33) 


~ 9 - 




































The bulk term vanishes by Eqs. (|2.28), (2.31), and the boundary term gives the on-shell 


action. We therefore have from (2.11) 


q{t,x)= lim +2A'^{t,x;u)) = lim [-KGtxit,x-,u)]j^^^^ , 

XL ^OO XJj ^OO 

(2.34) 

and 

nA{t,x) = Nc lim [^/T^F*“(^, x;«)], jA{t,x) = W Hm [y/^F^^{t,x]u)\^^^ . 

XL aOO XJj ^OO 

(2.35) 


In (2.35), we have introduced a “bulk axial current” [72]: 


Mu) = ■ 


(2.36) 


As in general bulk current J^(u) and bulk held Gtx{u) might be divergent near the bound¬ 
ary ti —)• oo. we use the subscript “Ren” in ( 2.35|) and ( 2.34 ) to denote the subtraction 


of such divergences in ( 2.34 ). The correspondence: ( 2.35 ) and ( 2.34 ) has been used in 
Ref. [52]. It is also interesting to note that from the u-component of the hrst equation in 


(IllSal) : 


Ncdfx = 2KGtx , 


(2.37) 


the anomaly relation ( [1.2| ) will be reproduced by taking tt —>■ oo limit on both side of ( 2.37| ) 
using (|2.35|) and (|2.34|) . 

We would like to comment on the nature of current Ja- Naively Ja obtained by a 
functional derivative is by dehnition a consistent current with respect to havor gauge. In 
fact, it is also the covariant current. We can conhrm this by noting that boundary source 
entering the bulk held strength only through boundary values of Ea and G, which are 
axial gauge invariant, therefore ja is manifestly axial gauge invariant. The agreement 
of consistent and covariant currents may appear odd: this is because the QCD anomaly 
studied in this section is realized with an on-shell action that is manifestly axial gauge 
invariant. Therefore, the current obtained from functional derivative is also invariant, as if 
it were an ordinary current. In contrast, the QED anomaly is realized with an anomalous 
on-shell action under axial gauge transform. In this case, we can not have a current which 
is both conserved (consistent) and invariant (covariant), in the presence of external axial 
held. As an example, we will see that the covariant current is not conserved in section. Q, 
where we expand our study to include QED anomaly. Eurthermore, because holography 
has access to gauge invariant (with respect to SU{Nc) gauge) quantities only, ja is also 
covariant current with respect to SU{Nc) gauge transform. We stress that the action ( 2.32| ) 


is very different from what we would have obtained by a naive substitution of ( 2.26 ) into 


( 2.22 ). In particular, the kinetic term of M would have an opposite sign, which would lead 
to a wrong sign for q, [48,69,70]. 

2.4 Computing one point function 

We now ready to compute one point function. We will work in Fourier space: dt —>■ —iw, 
dx —)• ik. Then using ( 2.37|) and Bianchi identity, one would express bulk current J\, 
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as: 


J\{u;uj,k) = 
J%{u;u},k) = 


Ncy/—^{ik)fE'j^{u] (jj, k) — 2K{iuj)G{u; uj, k) 
— k? f 

Ncy/—^{ioj)f Lo, k) — 2K{ik)fG{u] to, k) 


— k'^f 

Here we have introduced short-handed notations: 

Ea{u; uj, k) = -Ftxiu] uj, k) , G{u; uj, k) = -Gtx{u; uj, k) 


(2.38a) 

(2.38b) 

(2.39) 


Here and hereafter, we use prime to denote the derivative with respect to u. 

From (|2.38|) , we observe that we only need to solve equations for G{u; uj, k),EA{u; uj, k) 
to obtain one point function nA,jA,Q- From (2.25|), one finds: 




G"+ — + 


u /(a;2 - k^f) 


a + k^AkIla = 


v?p 



■ 4a:2 ■ 

G + 

2Kujkf' 

UcJ 


[p^f (a;2-A:2/)J 


Ea , (2.40a) 


E’U 


uJp' 




2u /(a;2 - kkf) 

1 

N, 


. (oj^ - k^f) ^ 

Ea 4 -— Ea — 




4^2 




Ea + 


2Kujkf 


E'cJ - ^V)_ 


G. (2.40b) 


It is understood that the back-reaction of the flavor branes will induce 1/Nc correction 
to the black-brane metric. Analysis shows that the correction to the metric could induce 
terms ~ G/N^ and ~ Ea/Nc to ( 2.40a ) and terms ~ Ea/N^ and ~ G/N^ to ( 2.40b|) . we 
will seek solutions to the leading nontrivial order in power series of l/Nc- 


G = + ... , 


EA = Ef + —E)l' + .... 


Ak 


Nc ' ^ ^ ' N, 


(2.41) 


To the leading nontrivial order, the solutions are not affected by the back-reaction. Ac¬ 
cordingly, we will compute one point function such as q,nA,jA to first non-trivial order in 
Nc- From ( |2.25| ) and ( p.41 ), we found that E'A'^ satisfies the homogeneous equation: 

uJp' 


E^A^" + ( :^ + 


^2u ' f{uj'^ — A:2 /)) ^ 

while G^^^ satisfies in-homogeneous equation: 






(2.42a) 


G«.)» + (_i + ff \ G(0)' + td^Gl") = , ,2.42b) 

V u fiuj^-k^pJ vAp Ggf{uj^-k^f) A ^ y 


At order 1/Nc, we further have: 


E^P" + 


Pf 




2u - A;2/) 


+ 


{P-kp)BA ^(1) 


-kpr^ ’ 

2Kkujf' 


EAp 


_ _ 

“ G^u^/p{uj^ -kp) 


(2.42c) 
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and a similar equation for 

Behavior of lo, k), E^^\u]U],k) and E^^\u]U!,k) near boundary can be deter¬ 

mined directly from (|2.42|): 


G^^\u; uj, k) = a 2 (w, k)v?{l -!-•••) + ^)(1 + ''') > (2.43a) 

Ef{u-,uj,k)=EQ{u:,k){l + ---)+Ef\u:,k)(u-^/^ + --^ (2.43b) 

E^);^ (n; w, k) = E^'^ (w, k) + • • •) , (2.43c) 

where ... denote terms higher order in 1/u. In ( 2.43cj ) we have defined the solution to 
in-homogeneous function ( p.42c ) e'^^ in such a way that E^^^ {u —>■ oo) = 0. 02 here is 
related to at, ax, 9 by 


02 = 


K 


[(w^ — k‘^)9 + 2iojat + 2ikax] 


(2.44) 


In deriving (p.44), we have used the relation: 


K 


u ^duG = [f ^dt {dtM 2At) - dx {dxM + 2Ax)] 
zu. 


(2.45) 


which can be derived from ( |2.3C| ) and (2.31). It is useful to note that £'o )®2 are invariant 
under transformation ( |2.15| ). 

As usual, we impose the infalling wave condition at the black hole horizon for (2.42|): 


lim G^^'^ {u] ui, k), E^?\u-, ui, k), E^}\u; cj, k) ^ {u — 1) . 


U^Ufi 


(2.46) 


Here uh = 1 denotes the location of horizon. Consequently, with given boundary value 
02 , Eo, ( 2.42|) can be solved and bo{uj, k),E^'^ (ca, k),E^^^ (ca, k) will be determined from the 
resulting solutions. They are related to one point function q, nA,jA via (2.38) and definition 
(|2.35D , (|2.34D . We therefore have 


q{uj,k) = Kbo{u],k ), 


nA{u],k) 


jA{^,k) 


1 

a;2 — k"^ 

I 

a;2 — 


-^{ik) {N,Ef\oj,k)+E^^\u,k)) 
{N,Ef\uj,k) + E^^\u,k)) 


2K{iuj)bo{u}, k) 
2K{ik)bQ{uj, k) 


(2.47a) 

(2.47b) 

(2.47c) 


3. Medium’s response to chiral charge imbalance 


In section, we will solve ( ^.42 ). with two different boundary conditions and consider the 
relation between Ja, q and ua- Physically, we would like to use those two different boundary 
conditions to model two different mechanisms for the generation of axial charge imbalance. 
In particular, we consider: 
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Case 1 axial charge imbalance is generated by a domain with non-zero winding 
number. To model this situation, we set axial gauge field to be zero at boundary, 
i.e., at, ax = 0 but turn on a non-zero 0{u}, k). Consequently, boundary condition for 
(|2.42|) becomes: 


02(w, k) 



k^)9{uj,k) , 


Eq{ui, A;) = 0 . 


(3.1) 


Case 2 axial charge imbalance is generated by non-topological fluctuations. To 

model this situation, we instead set 02 = 0 and consider a non-zero axial electric field 
on the boundary: 

a2(w,A;) = 0, Eo{uj,k)j^O. (3.2) 


3.1 Medium’s response to axial charge imbalance generated by topological fluc¬ 
tuations 


In this section, we will study medium’s response to axial charge imbalance generated by 
topological fluctuations. As we discussed previously, this amounts to solve ( 2.42|) with 
boundary condition ( ^.l| ) and ( 2.43cD . As in this case, there is no source term for ( ^.42a ), 
E^^\u‘,u!,k) = 0 trivially satisfies ( |2.42a| ) and consequently ( 2.42b ) becomes an homoge¬ 
neous equation: 

We will seek the in-falling solution for (3.3) in hydrodynamic regime w, A; <C 1. In this 
regime, the solution can be obtained analytically by first solving (^) order by order in 
power of w, A; away from horizon and then determining integration constants by matching 
with in-falling wave boundary conditions near the hrozion. Away from the horizon, we can 
drop the third term in (3.3) to obtain the following solution 


lUJ 


gh{u) = 1 - (-^ 


where we have defined 


ru 

/ ' 

JUH 


du' 


3Cg(l - sV) 

k 


1 


u' — 1 


+ log(u - 1) 


s = 


u 


to save notations. It is easy to check that behavior of (^) near the horizon u 
matched to the infalling wave behavior in small uj limit: 


(3.4) 

(3.5) 
1 can be 


(u — 1) *3=1 — — log(u — 1) -I- O(w^). 
3 


(3.6) 


It is also worthy mentioning that the integral over u' in (3.4) is convergent as we have 
explicitly taken the log(u — 1) outside the integral. 

From boundary condition (p.l|) , we then fix the normalization of : 


G^^\u;u},k) = C-^e)gh{u-,u},k). 


(3.7) 
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Expanding (^) near the boundary, we obtain 
q{uj,k) = 2Kbo{uj,k) = 


/A'2^ 


\ r , Ca n^l 


{iuj — k lim 

Ei—>-00 


\Gg ) 

JuH '\/~9 2 


+ 0{u}^)}9{u},k). 


(3.8) 


In (3.S), the subtraction is necessary to remove the divergence near the boundary. As 
the ratio —q{uj, k)/9{u}, k) should be matched to the behavior of retarded Green’s function 
(0), we identify Ecs, in the present model: 


Tcs 2iE2 2iE2^3 


C„ 


Cn 


rp6 _ 


8X^T^ 

7297rM|,^ ’ ’ 


Kcs = -2 


(lim [ r du'u' 

\Cg) VJuH 2 J 2 V 


bcs 






2437r2 


(3.9) 

(3.10) 


where in the last step, we recover the units and used (p.24| ) and definition ( |2.1^ ). The Tcs 
in Sakai-Sugimoto model was computed previously in [73]^ (see Ref. [74] for Tcs in other 
holographic models). 

To compute nA,jA, one needs to solve in-homogenous equation (2.42c): 


+ ( X- + 




2u - k^f) 


E<i>' + = 

A A 


2Kkujf' 


— k'^f) 


(3.11) 


with G^^\u;ui,k) given by (^). At leading order in uj,k, the solution reads 

EP{u-u;,k) = -(^^^Cg^9{co,k)u-^/pl + Oico,k)) , (3.12) 

which can be easily verified by substituting it into ( 2.42cD and comparing results at leading 
order in io,k. Now substituting (|3.7[) and (^.12 ) into (|2.47|) , we have obtained the axial 
charge density generated by 9, 


UA = 


2K^ 

“a 


9{u},k)[l + 0{u},k)] = 


Ecs 


9{uj, k) 


(3.13) 


and jA vanishes at this order. To obtain j^, we consider t,x components of ( 2.25b| ) 
in the presence of given by (3.7). We then obtain flow equations of bulk current 
Ja{u; uj, k), TKu; cj, k) along radial direction u: 


duJi = -2KG^P + {ik^g^^g^^) = 2KGP + {ioj^g^^g^^) E 


(0) 


n(l) 

'A ■ 

(3.14) 

Now using the relation between G^J and GpP G^ which can be obtained from ( 2.25aD in 
the absence of E, 

AO) 

(3.15) 


^( 0 ) _ -ikfd^Gl^J 

'-'ut — ' 


UJ‘ 


-ky ’ 


(0) ^ i^duG^ 

>-7 


UJ‘ 


-k^f 


^our results (3.9) has a different normalization from [73] 


- 14 - 













































we 


have from (^): 

d^j\ = + {ikV^g^'g^^) 4'^ = [1 + Oiu^, k)] . (3.16a) 

2ikK^6 , . - ii fi) 2ikK‘^6 , / . x 

—^=- + (za;v^5 ) 4 ^ =- 1=- [1 + 0{uj, k)] , (3.16b) 


duJ% = - 




V^g 


By integrating over u, we therefore have: 

nA{uj,k) = nA{uj,k) + AnA{uj,k), jA{uj,k) = jAiuj,k) + AjA{uj,k). (3.17) 

Here n^,j^ are values of bulk current J\{u] ui, k), J^{u; uj, k) at horizon u = uh- We 
already know ua 0(1) from (|3.13|) and 


Aua = {—2iuK^6) 
Therefore we must have 


- 1-00 2 

.4 


0{uj,k). 


- 

Ua — 


2K^ 


9{uj,k) [l + 0(a;,A:)] 


On the other hand, 
AJa = i-2ikK^e) 


/■OO . 

/ 

Juh 


Ren 


-2ikK‘^e, , 

= (— n—> 4^ 

tyg U^OO 


LJ Uff 'J 5 _ 

By comparing (|3.2C|) with ( p.lO|) , we have 

Aja = Kcs{ik)0{u,k) 


■ pOO 

/ du'u' — f 2 
Juh 


(3.18) 


(3.19) 


(3.20) 


(3.21) 


The current on the horizon can be determined by substituting ( ^.12 ) into ( 2.38 ) and 
taking u —?> uh limit: 


Ja = lim [>/l(w;w,fc)] = lim 

U^Uh U—^Uh 


y^{iuj)fduE^j\u-,uj,k) - 2K{ik)fG^^\u-,uj,k) 


up' — k'^f 


= lim 

U—^Uh - 


c\ T^Z 

V—74^4; k) =—^r^{ik6). (3.22) 


30„ 


Here we have used a property of any function satisfying in-falling wave boundary condition, 
say Zin(u) that 

lim (fduZiP) = —iuj lim Zm . (3.23) 

U—yUfj 

Comparing ( 3.1S| ) and ( 3.22 ), we found that on the horizon, Ja and are related by 
Tick’s law: 

if = -OVnf . (3.24) 


To establish (3.24), we also used the value of diffusive constant D (3.30) in current model. 
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To sum up, in this subsection, we have studied axial current in response to axial 
charge imbalance created by topological fluctuations. To represent a domain with non¬ 
zero winding number, we first turn on a non-zero k) and found that it will induce 
a non-zero q ( |3.8D and consequently a non-zero axial charge density nA{(-o,k). The axial 
charge density nA{uj,k) and 6{co,k) are related by ( 3.13| ). Furthermore, the induced axial 
current can be divided into two part. The first part is due to the diffusion of ua while the 
second part is in the opposite direction to the diffusive current and is proportional to kqs , 
which quantifies the kinetic energy carried by a topological domain. We verified relation 
B as first proposed by us in Ref. [52]. It is interesting to note that holographically, 
the diffusive (dissipative) current coincides with the current on the horizon (^.24) while 


the non-dissipative current (O) is given by the integration from horizon to the boundary 
(c.f. (^3oD and (|T2ID ). 


3.2 Medium’s response to axial charge fluctuations generated by non-topological 
fluctuations 

We now consider medium’s response to axial charge imbalance generated by non-topological 
fluctuations. Following our discussion in Sec. 2.4 we will solve ( p.42 ) with boundary condi¬ 
tion (3.1) and (2.43c). We first need to solve the homogeneous solution (2.42a). Similarly 
to (^), the infalling wave solution to ( 2.42aD reads: 

,f3C^{l-s^f) 1 


eh{u-,uj,k) = 1 - 


lOJ 


du 


uh 


V^f 


u' — I 


+ log(u - 1) 


Consequently from boundary condition ( |3.2D , we have: 

'Eo{uj,k) 


E^^\u-,uj,k) = 


eh{u]UJ,k) , 


_CEi0J,k) 

where ce{oj-, k) is defined by the value of eh{u;u, k) at boundary: 


(3.25) 


(3.26) 


ce{oj, k) = eh{u ^ oo; u,k) = 1 + iio 


OO ^ 

2 / W 


J Ul 


s~ I du 

UH 


-7 


+ 0 ( 1 ) 


= l + ^{s^ + 0(l)) , 
(3.27) 


Plug (|3.26| ) into ( |2.47]) , we obtain: 

• c - AT n ^o(w, A:) 
jAy^i k) — NqC^ 


nA(uj, k) = NcS Cr 


Eoiui^k) 


(3.28) 


CE{U},k) ’ CE[0J,k) 

The conductivity cr, diffusive constant D and susceptibility y can be extracted from ( 3.2j) 


as follows. First of all, in the homogeneous limit s —>• 0 of (3.28), we will reproduce Ohm’s 
law jA = ctEa- Therefore: 

u - ■ (3-29) 

On the other hand, Eqs. ( tl.28 ) must have a hydrodynamic pole corresponding to diffusive 


mode at w = —iDk^. This implies that ce{ui = —iDk'^, k) = 0 hence: 

/•OO /O 

./ ^7 


D= du 
JuH 


-7 


r du'{u')-^/^=\f-^= ^ 

J Uh 


2txT 


(3.30) 
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In (^) and ( pO|) , we have used expression ( ^. 281 ) and have recovered the units at the 
last step. Finally, using the Einstein relation a = x^, we obtain the expression for x from 
(|339D and ( 1 ^): 

n T roo 

(3.31) 


X"' = - = — 
cr Nr 


f 

J Uf 


du' 


-c JuH V 7 

This relation between x and the bulk integration over -^—7 is in agreement with general 
expression in Ref. [72]. 


Eq. ( 3.28 ) implies that turning on an external axial electric field Eq would generate 
a local axial density and axial current. It would also create a non-zero q, which can be 
determined by solving (p.42b): 


(7(0)// + ( _ 


with given by (|3.26). At leading order in the in-homogeneous solution reads: 


\ ^(0)/ ^ 
^/(a;2-fe2/)y 


2Kku}uf ( 0 ) 
Cgf{io^-kp) ^ ’ 


(3.32) 


G^°\u;u,k) = 


w (1 -|- O{oj, k)) H—: 


2gh{u-,u,k) 


iu!{l — s2) 


Eo{uj,k) 
CE{ui,k) ' 


(3.33) 


As one can check, the first term, i.e. , sKv?/Cg term is a special solution to in-homogeneous 


equation ( 3.32 ) at leading order in w, k. gh{v) (|3.4|) , the solution to homogeneous equation, 
is introduced to guarantee boundary condition p. 2 |). As a result, we have: 


q{u,k) = 


f2K^\ 

—is 

EQ{uj,k) 

f2K^\ 

—is 

EQ{uj,k) 

V c, ) 

(1 — s2)a; 

CEiuJ,k) 

V Cg ) 

—s'^uj 

CE{uJ,k) 


(3.34) 


In the last step, we dropped the 1 in the bracket. This is justified in the diffusive regime 
where P ~ u. To hnd the response of q to ua, we first note that q in (3.34) contains 


responses to both ua and axial electric held Eq. In fact, we can exclude the response to the 
latter by considering ua induced by a normalizable mode. This occurs when ius^ = —3/2. 
According to ( |3.27 ) and (3.28), it implies that ua remains hnite while both ce and Eq 
vanish. Physically ua in this case is induced by a diffusion wave. The response of q to ua 
is then given by 

q 2K^ 1 4^2 

UA 


NcCgCj — ZWs2 


3NcCgC.y ’ 


(3.35) 


where in the last step we used = —3/2 and dropped higher order terms in uj. We note 


that (|3.35|) is precisely (|1.5|) as we advocated in the introduction. 

We now show that the intuitive argument above could be established more rigourously 
in real space. Eor this purpose, it is convenient to perform the Eourier transform over uj 


and directly consider nA{t,k),jA{t,k) and q{t,k). By dehnition and ( 3.28 ), we have: 
jA{t,k) 


Nr. 


f 


duj 


e *‘^7 'a(w. A:) = 


= -a 


f 


r 

ujEq{uj, k) 

1-00 27r 

UJ -|- iDk"^ 


duj 

—( 
27r 


—iujt 


UJ 


UJ + iDk"^ 


dt'E^^'EQ{t\k). (3.36) 
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The above integral is nonvanishing for t > t' when we can pick up the diffusive pole in the 
lower half plane. We then obtain 

= -C^Dk^Ioit, k). (3.37) 

/o(A:) here is dehned by: 

Io{k, t) = = -Wk^, k)= j k). (3.38) 

J —OO 


A similar computation gives: 


nA{t,k) 


Nr. 


= {-iC^)kIo{t,k), 


,2K^ 


q(t, k) = {-iDk) Io{t, k). 

0/7 


(3.39) 


(3.40) 


In (|0^) , we have neglected a contribution higher order in k. 

We assume the external field Eq only exists in a hnite time window. At sufficiently 
late time t, we can regard ja and q as responses to ua- Comparing ( 3.31: ) and ( 3.39|) and 
return to real space, we obtain Tick’s law: 


jA{t,x) = -DVnA{t,x). 
Moreover, ( 3.391) , (|3.40|) lead to the relation: 


(3.41) 


q{t,x) = 


2K‘^ 


D 

c2 


nAit,x) = 


res 


nA{t, x) = riA{t, x) = , 

(3.42) 


where we have used (P), P). 

The results of this section can be summarized by the following response matrix, char¬ 
acterizing the response of q, ua and Ja to 9 and Ea- 


( q ^ 

( 

UA 

{Ja) 

- 


ik 


ituT'es 

2T 

Tes 

T 

I^CS - 


T k(Lj^iDk‘^) ' 
ka 


ua 

u+iDk'^ 


/ 



(3.43) 


We keep only terms lowest order in small u) and k limit. We further restrict ourselves to 
the regime co ^ k for the case with source 9, and to diffusive regime ~ cu for the case 
with source Eq. The transport coefficients we presented are their first nonvanishing order 
in Nr', the responses of ua and Ja to Ea are at order 0{Nc), while the rest responses are 
at order 0(1). 

Before closing this section, we would like to comment on the 0(1) correction to the 
responses of ua and Ja to Eq above. This requires us to go beyond leading order in 1/Nc 
and compute E^^'^ from ( p.40b ) with E^'^ given by (3.2f:). Similar analysis shows that near 
the boundary, E^^'^ ~ + • • • as u —)■ oo, which would give divergent contributions to 
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UA and Ja- This is because the mixing of the bulk fields changes the dimension of the 
operator. We note that the change of operator dimension occurs immediately with mixing 
in bottom-up model in [71], while in our case it occurs from the subleading order in l/Nc- 
As we explained earlier that the solution this order in 1/Nc is incomplete without including 
back-reaction of the flavor branes. It is curious to see if including such backreaction would 
remove the potential divergence. Although these higher order corrections do not affect the 
results of our paper, we hope that we could revisit the puzzle in future. 


4. Chiral Magnetic Effect and universality 


In the previous section, we have considered two different situations where axial charge 
imbalance is generated, we now want to study whether the CME current would depend 
on the microscopic origin of axial charge density. In particular, we will compute the ratio 
between CME current and axial charge density ua in low frequency and momentum 

limit as defined in (lA). If axial charge density is static and homogeneous and CME current 
is universally given by (M), one would have: 


Xdyn — X 1 


(4.1) 


due to linearized equation of state Sua = However, it is not obvious if (4.1) would 

still hold if UA is generated dynamically as considered in this paper. Of particular interest 


is the case considered in Sec. ^ that axial charge imbalance is created by topological 
fluctuations. 

To compute O), we turn on a small background magnetic field Fy^ = eB, i.e. mag¬ 
netic field is longitudinal to the direction of in-homogeneity as considered in the previous 
section. Then in the presence of bulk axial field F^, F^^ and F^, F^, F^^,F^ and F^ 
components of vector field strength will be excited due to Wess-Zumino term: 


Awz = Cs A tre-^«/ 2 ^ - C3 A , 


(4.2) 


the same form as that of axial gauge field A and WZ term from (4.2): 


with F 4 = dC^ given in (M). The action for vector field consists of DBI term, which has 

(4.3) 


Sy = -Nc J d^xdU 




where 

Here, we use subscript/superscript V for vector gauge field strength (below we will 
also use subscript/superscript A for axial gauge field). 

The variation of Sy gives the equation of motion: 


Om = Kse^^^F^Q . 


(4.5) 
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^MQ taken from solutions obtained in previous solutions. To compute (IT), it is 

sufficient to work at linear order in eB. We therefore could neglect back-reaction due to 
eB to the holographic background and solution obtained in the previous section. 

As before, we define the bulk vector current 


J^{t, x; u) = {t, x; u) 


(4.6) 


One point functions nvit,x),jv{t,x) are similarly given by the boundary values of bulk 
current: 


nv{t,x)= lim Jy{t,x;u), jv{t,x)= lim Jy{t,x;u). 


(4.7) 


As before, the vector current defined here is a covariant current. Similar to (2.36), it would 
be convenient to express vector current in terms of Ey and Ea as 


Mu) = w 


{ik)y^fE'y{u) - {2i(jjKB)EA{u) 


J^{u) = N, 


— k‘^f[u) ’ 

{-iu)y/^f E'y{u) - {2ikKB)f Ea{u) 


a;2 — k‘^f[u) ’ 

where we have introduced the short-handed notation for “bulk electric field” 


Ey{u]UJ,k) = -E^{u;uj,k). 


(4.8a) 

(4.8b) 

(4.9) 


We can easily verify using (O) that the covariant current is not conserved in the presence 
of external axial field Ea- The equation for Ey{u;u},k) reads 




E!^+{— + 


2u 

2Kb 


a;2 - fe2/ 


_ M - k'^f) ^ 

H 7372 - 


[u®/2/ (^2 — k^f) 


U^f 

+ Tr 


2KLokf 




G} + 0{{eB)^). (4.10) 


We will solve ( [4.100 with Ea,G determined in the previous section. To concentrate on 
vector current induced by axial charge imbalance, we will not turn on any source for vector 
field, i.e. , imposing Ey{u -A oo) = 0 on the boundary and use the standard in-falling 
wave boundary condition on the horizon. 

To compute the vector current, it is also convenient to write down “flow equation” for 
bulk vector current by taking t,x components of ( |4.5| ) and using definition (4.6): 

duJviu) = ^^ Ja{u) - Ey{u) = ^^ Ja{u) [1 + 0{uj, k)] . (4.11a) 


2Kb 


yMlf 
2Kb rt 


duJviu) = -j=Ja{u) - {iujy/^g^^g'^^) Ey{u) = --j=Ja{u) [1 -F 0{uj, k)] , (4.11b) 

On the R.H.S of (gUbj) , we have used the fact that Jy^ term is always dominated over 
Ey term in small oj,k limit due to additional gradients in front of Ey. This is because 
from ( 4.10 ) we observe that Ey is the same order as Ea and from ( 2.36 ), is at least 
the same order as Ea- 
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Now integrating ( 4.11a ) over u and using definition (^), we have: 
nv{uj,k) = Anv{u},k) + ny{uj,k), jviuj,k) = Ajv{uj,k) + jy {uj,k). (4.12) 


Here Uy and jy are values of bulk current Jy{u;u,k) and J^{u;uj,k) at horizon u = uh 
and 


Any{uj, k) = Kb 


■ poc 

-JuH 




Ajv{u}, k) = 2Kb 


/ 




.13) 


We now claim that CME current should be identified with Ajy, i.e. , 

,j\{u-,uj,k)' 


jCME ^ 


/ 


du 


\_^UH 


-7 


(4.14) 


The physical motivation behind identification (4.14) is that generically in a holographic set¬ 
up, the current on the horizon is dissipative (see also example below). On the other hand, 
the CME current is non-dissipative. Therefore one should exclude the horizon current from 
the total current when identifying CME current holographically. 

We now consider the implication of (4.14). With (f4.14) and (p.31), Xdyn becomes: 


Xdyn = lim {nA(w,A:)/ 
uj,k^0 


f 




} 


It is clear that if in small ca, k limit, bulk axial current is constant, i.e., 

j\{u] ui, k) = nA{oJ, k) [1 + 0{ijj, k)] , 


(4.15) 


(4.16) 


it follows one can replace J\{u]OJ,k) with ua in ( 4.15| ). Consequently, one will arrive at 
(4.1) by noting ( ^.311) . Therefore ( [4.16| ) can be interpreted as a condition for the validity 
of (p|). 

In both cases considered in Sec. |^, the condition ( 4.16 ) is indeed satisfied, we therefor 
have (0) for those cases. 

Eor completeness, we will calculate total jy for both cases. Eor the hrst case (c.f. Sec. 3.1), 
it is straightforward to check that jy ~ 0{uj,k)6, which is sub-leading compared with 
we therefore have: 


jy{uj,k) =j\j^^{uj,k) = ( — 

X 


re 

bcs 






(4.17) 


where we have used (0)- By comparing ([4.17| ) with (|1 . 1|) , it is tempting to make the 
identification: 


fkA = 


2rsph 


(4.18) 


In Ref [4], ^a is identified with dtO. ( 4.18 ) corresponds to replace dt in dfO with the inverse 
of characteristic time scale of sphaleron transition l/rgph. 


- 21 - 














































The situation is different for the second case (c.f. Sec. ^). In this case, jy is the 
same order as jy^^ in small u, k limit. From (|3.28D and (4.15D, we have: 




2NcKb \ f Eo{uj,k)\ 4 :NcKb { kEo{uj, k) 


X 


CE{u),k) 


\LJCe{!^, k) 


(4.19) 


On the other hand, to obtain nv,jv, we need to solve (|4.10D , at leading order in Nc, that 

2KBUjkf' 


E'{r + 


5 uj^f 
+ 


p/ I 


2u ' u)‘^ — k?f) ^ ^ f {uj‘^ — k^ f) 

where E^ is given by (3.26[). The leading order solution to (|4.2C|) reads: 


Ev = - 


iKBEo{u;,k) 
3ce{ui, k)C^ 


i-3/2 + o{u}, k) 


s. 


Ea, (4.20) 


(4.21) 


Again, ( |4.21 ) can be easily verified by direct substitution. Now substituting ( [4.21 ) into 
), we have: 

Eo{uj,k) 


nv = —2iNcKB- , ,. , 

ujce[co, k) 

and jv vanishes at this order. Since Any ~ E){1), we thus have 

Eo{uj,k) 


ny{u), k) = nv{uj, k) = —2iNcKB 


u)Ce{oj, k) 


(4.22) 


(4.23) 


jy is obtained in the same way as in the previous section: 


H ■ 


jy{u},k)= lim [Jy{u;u},k)] = lim Nc 

U—yUff 


{-iuj)yN:jfduEv{u] w, k) - 2{ikKB)fEA{u] u, k) 


U}‘ 


-k^f 


4NcKBkEo{uj,k) 


= lim —Nc \y/—^Eviu]io,k)] = - - —--—=—iDkny{co,k). (4.24) 

'll _i-l/ 7-T L JO/ - I, < 1yi\ 


U^Uh 


ujce{u}, k) 


Again we see that jy can be interpreted as a diffusive current. As in this example, both 
diffusive current jy (4.24) and CME current (4.1£) would contribute to the total vector 
current in small uj^k limit, to compute CME coefficient hence Xdyn properly, it is crucial 
to identify CME contribution, i.e., ( 4.14 ). 

To close this section, we would like to comment that while in this paper, we are working 
in a specific holographic model, the relation ( 4.12| ) still holds for holographic action for bulk 
vector field of the form (|4.3|). Consequently, assuming the identification of CME current 
( 4.14 ), the condition ([1.16|) would warrant that Xdyn = X- Moreover, the violation of 
condition (4.16) would also break the relation (4.1). 


5. Stochastic hydrodynamic equations for axial charge density 

We now formulate a hydrodynamic theory for axial charge density by including stochastic 
noise from both topological fluctuations and thermal fluctuations. We will focus on the 
dynamics of axial charge thus setting temperature and fluid velocity to be homogeneous 
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and time-independent. We could then work in the frame that the fluid is at rest: = 

(1, 0,0, 0). To close anomaly: 

dtnA{t, x)+V ■ jA{t, x) = -2q{t, x). (5.1) 

we want to express q and ja in terms of noise and ua (or its gradients). The constitute 
relation, which relates axial current ja to n^, is of the conventional form: 


jA{t, x) = -DVnA{t, x) + ^(t, x ), (5.2) 

where ^{t,x) encodes axial charge generated by thermal fluctuations: 

{^{t, £C)) = 0 , x)^j{t, x')) = 2aT6ij6{t - t')5^{x - x'). (5.3) 

Here (...) denotes the average over noise and i^j = 1, 2,3 run over spatial coordinates. The 
magnitude of ^(t, x) is given by the standard fluctuation-dissipation relation. Furthermore, 
q{t,x) can be related to nA{t,x) using (|L^): 

q{t,x) = -gg^(t,a;). (5.4) 

^q(t,x) is the noise due to topological fluctuations: 

{U(t,x)) = 0, {^q{t,x)^q{t,x')) = TcsS{t- x'), (5.5) 

and we will assume that there is no cross correlation between two different types of fluctu¬ 
ations: 

{^q{t,x)Ci{t,x')) =0. (5.6) 

This completely specifies our stochastic anomalous hydrodynamic equations. The noise 
due to topological fluctuations ( |5.5|) have been considered previously in Refs. [75]. On 
the other hand, for a conserved current (such as vector jy), the noise of the form ( p.3| ) is 
standard. Incorporating both fluctuations in ( ^.Sj) is new to the extent of our knowledge. 
As an application, we will consider equal time axial charge correlation function: 

Cnn{t, x) = {[nA{t, x) - n^(0, x)] [nA{t, x) - n^(0, a;)]) . (5.7) 


We start with equation for ua readily followed from the stochastic hydrodynamic equations 

dO): 


dt - + T, 


-1 

sph 


nA{t, x) = x) + 2(q{t, x). 


(5.8) 


Under the initial condition nAiO,x) = 0, Cnnit,x) characterize the magnitude of axial 
charge fluctuations at time t and location x due to (both) fluctuations. 

To compute (5.7), we first Fourier transform (|5.8|) into k space (but keep t-dependence): 


dt + Dk^ + r. 


Sph 


nA{t, k) = [-ik ■ ^(t, k)] + 2^q{t, k) 


(5.9) 
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The solution of under initial condition n/i(0,a;) =0 reads: 

r ^ 

nyi(t, k)= dt'e~^^^ +'^sph)h-*) [(-ifc) • k) + 2^g(t, fc)] . (5.10) 

Jo 

We therefore have: 

{nA{t,k)nA{t,k')) = f dti 

Jo Jo 

^ [hk'j{ii{tl,k)ij{t2,k')) + A{^q{tl,k)^q{t2,k'))] . (5.11) 

^) in Fourier space, 


Using (|5.5|) and I 

{i{t, fc)) = 0, (e*(t, k)iq{t, k')) = 2aT6,j6{t - t')6^{k + k '), 

{iq{t, k))=0, {Cq{t, k)Cqit, k')) = TcsS{t - t')6^ {k - k') , 

and performing the average over noise, we have: 

{nA{t, k)nA{t, k')) = 2 {aTk^ + 2rcs) / ^ 

Jo 

aTk'^ + 2rcs 


(5.12a) 

(5.12b) 


+ -^sph 


1 _ g“2(i5A:2+T^pl)i 


S^{k + k') = xT 


I — g“2(^fc^+T^ph)t 


5^{k + k'). 

(5.13) 


In the last step, we have used Einstein relation a = xD ^-nd (1^). Now returning to real 
space, we have: 


Cnn{t,x) = xT 


d^k 

w? 


Axk 


1 _ g“2(i5fc^+'r„p),)t 


= {XT) 


5^{x) - 


1 


{SnDt) 


3/2 


- _i£r 

g ■^sph g SDt 


(5.14) 

It is worthy noting that as we are in hydrodynamic regime, nA{x) here should be understood 
as the coarse-grained axial charge density inside a fluid cell and x is the spatial coordinates 
labeling the corresponding fluid cell. 

We now discuss the implication of ( 5.14 ). At very early time that Dt <C (therefore 
t <C Tgph) where Lceii is the size of a fluid cell, the Gaussian appearing in ( 5.14 ) essentially 
becomes a delta function and we then have: 


Cnn{t,x) K, xT 


2t 

1 — g 


^^{x) PS ^ 16(x) = dFcs i d^(x) 

Tsph 


(5.15) 


where in the last step we have used (1.6). At this stage, there is no correlation among axial 
charge in each fluid cell (c.f. the delta function in (5.15)). Integrating (5.15) over volume 
f d^x, we further recover relation between the fluctuation of axial charge and Chern-Simon 
diffusive constant: 

(Qi) =4rcsW, (5.16) 

where V is the volume of the system. While it has been widely used in literature to estimate 
the fluctuation of axial charge, ( |5.16 ) is no longer valid at the stage that Dt Lin- In 
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this stage, the spatial dependence of axial charge fluctuations in ( 5.14 ) become important. 
The diffusion generates additional spatial correlation among axial charge density. Finally, 


in the long time limit t S> Tgph, the second term in (5.14) are suppressed exponentially and 
axial charge fluctuations are given by: 


Cnn{t 00,x) ^ (xT) 5^{x) . 


and 


lim (Ql) 

E —>-00 


XTV; 


(5.17) 


(5.18) 


This is of course expected as in long time limit, the (Q 5 ) should approach its thermal 
equilibrium values (|5.18| ). 

Finally, we remark that to obtain ( 5.1^ ), we have used the relation between Chern- 
Simon diffusive constant Fcs and sphaleron damping rate Tgph O)- Therefore like Einstein 
relation a = xD connecting conductivity a and diffusive constant D, the relation between 
Fcs and Tgph is also fixed by the requirement based on thermodynamics (p.l8|) . It is reas¬ 
suring that the relation ( p.. 6 [ ) is also realized, as we discussed in Sec. 3.2, in the holographic 
model studied in this work. 


6. Summary and Outlook 

We have analyzed the anomalous transport of a non-Abelian plasma in a de-confined phase 
with dynamically generated axial charge using a top-down holographic model. In partic¬ 
ular, we consider two separate cases in which the axial charge is generated due to a) 
topological b) non-topological thermal fluctuations. When the axial charge is generated by 


topological gluonic fluctuations, we show a non-dissipative current (1.3) is induced due to 


chiral anomaly in Sec. 3.1. We also illustrate holographically the damping of axial charge 
due to the interplay between flavor sector and gluonic sector. Furthermore, we consider 
the ratio of the CME current to the axial charge density at small ui and k (c.f Ql.?] )). 
We interpret such ratio as (the inverse of) “dynamical axial susceptibility” Xdyn (c.f. the 
discussion below (|l.7l )). We found in the context of current holographic model, dynamical 
susceptibility Xdyn is independent of the microscopic origin of the axial charge and coincide 
with the static susceptibility y. One phenomenological implications of our work, in partic¬ 
ular, Sec. Q is that axial charge generated by topological fluctuations and non-topological 
thermal fluctuations would contribute to CME signature in heavy-ion collisions. Eor this 
reason, we propose a stochastic hydrodynamic equation of the axial charge where we incor¬ 
porate noise both from both fluctuations in Sec. We found that the magnitude of axial 
charge fluctuations depend on the time scale that such fluctuations are measured (c.f ( 5.11 ), 
( ]5.18 )) and as (^.14 ) indicates, the diffusive mode would induce spatial correlations among 
axial charges. 

There are several issues that can be further studied based on the above analysis. The 
axial charge response to gluonic fluctuations can be studied when the flavor degrees of 
freedom back-react to the glue part of the theory. In this case, the coupling of the axial 
current and the gluonic topological operator, q, is not suppressed. 
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In this work, we found in hydrodynamic limit, the “dynamical axial susceptibility” Xdyn 
dO) is universal. It is both theoretically interesting and phenomenologically important to 
extend the definition of Xdyn and study its independence on finite uj, k. In this case, it is 
possible that the resulting Xdyn would depend on the origin of axial charge imbalance. 

In computing axial charge density correlation function Cnn ( |5.7|) from stochastic hy¬ 
drodynamic equation formulated in Sec. we consider a system in the absence of magnetic 
field. Once there is magnetic field, axial charge would also be transported by chiral mag¬ 
netic wave [51]. Furthermore, a new diffusive model would emerge due to the interplay 
between chiral magnetic wave and sphaleron damping [71,76] . It is interesting to see how 
those new modes would contribute to correlation among axial charge densities within the 
framework of stochastic hydrodynamics. 
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